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ABSTRACT 



We have used the high sensitivity and resolution of the IRAM interferometer 
to produce sub-arcsecond ^^CO J— 2-1 images of 9 protoplanetary disks surround- 
ing TTauri stars in the Taurus- Auriga cloud (7 singles and 2 binaries). The im- 
ages demonstrate the disks are in Keplerian rotation around their central stars. 
Using the least square fit method described in Guilloteau & Dutrey (1998), we 
derive the disk properties, in particular its inclination angle and rotation veloc- 
ity, hence the dynamical mass. Since the disk mass is usually small, this is a 
direct measurement of the stellar mass. Typically, we reach an internal precision 
of 10% in the determinations of stellar mass. The over-all accuracy is limited 
by the uncertainty in the distance to a specific star. In a distance independent 
way, we compare the derived masses with theoretical tracks of pre-main-sequence 
evolution. Combined with the mean distance to the Taurus region (140 pc), for 
stars with mass close to 1 M0, our results tend to favor the tracks with cooler 
photospheres (higher masses for a given spectral type) . We find that in UZ Tau E 
the disk and the spectroscopic binary orbit appear to have different inclinations. 

Subject headings: Stars: T Tauri - circumstellar matter - pre-main sequence - 
individual: BP Tau, CY Tau, DL Tau, DM Tau, GG Tau, GM Aur, LkCalS, 
MWC 480, UZ Tau - binaries: close - Radio-hnes: stars - 



1. Introduction 

Nearly all of our knowledge about the masses and ages of low mass young stars 
comes from their location in the HR diagram relative to theoretical calculations of stellar 
evolution to the main sequence. Despite considerable advances over the past 5-10 years in 
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understanding the structure and atmospheres of stars of mass M< 1 Mq, comparison of the 
currently available predicted evolutionary paths of young stars shows obvious differences. 
Empirical tests of the calculations have not been possible until recently because astronomers 
have not had independent measurements of either the mass or age of a young star. This 
situation is changing rapidly. One way to test the calculated tracks is to investigate whether 
they yield the same ages for stars expected to be coeval on physical grounds. Hartigan 
et al. (1994), Casey et al. (1998), and White et al. (1999) have apphed this test to wide 
binaries, the TY CrA system, and the GG Tau system, respectively, using tracks available 
at the time. Over a mass range extending down to the brown dwarfs, they found significant 
differences in the extent to which the several theoretical calculations satisfy the coevality 
requirement. 

The capability of mm-wave interferometers to resolve the spectral line emission of the 
outer disks of PMS stars offers the possibility to map their rotation. Since the disk mass 
is usually very small compared to the mass of the star, this provides the means to measure 
the stellar mass dynamically (e.g. Dutrey et al. 1994, DGS94; Guilloteau and Dutrey 1998, 
GD98). We use this technique to obtain new measurements of the mass of 5 single PMS 
stars and one PMS binary. Combined with our previous measurements of the PMS single 
CM Aur and the binary GG Tau (Dutrey et al. 1998, Guilloteau, Dutrey and Simon 1999) 
the measured stellar masses span the range ~ 2 to 0.5 Mq. We use the results to test the 
theoretical calculations of PMS evolution. Sect. 2 describes the sample of stars and the 
interferometric observations. Sect. 3 summarizes the analysis and presents the measured 
masses and related parameters. Sect. 4 compares the measured masses with those implied 
by location of the stars in the HR diagram relative to the tracks calculated by D'Antona 
and MazziteUi (1997, DM97), Baraffe et al. (1998, BCAH), Palla and Stabler (1999, PS99), 
and Siess et al. (2000, SDF). 
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2. Sample of Stsirs and Interferometer Observations 

2.1. Sample of Steirs 

All the stars studied are in Taurus- Auriga (Table 1). We obtained new interferometric 
observations in the ^^CO J= 2-1 line of the singles CY Tau, DL Tau, DM Tau, LkCalS, 
and MWC480 (HD 31648), and the spectroscopic binary UZ Tau E (Mathieu et al. 1996). 
They were selected because earlier interferometric observations (Table 1, Col. 5) indicated 
that they were, or were hkely to be, associated with resolvable circumstellar disks. BP Tau 
was observed because its HIPPARCOS distance, 56±14 pc (Favata et al. 1998), is very 
different from the 140 pc average distance to the Taurus star forming region (Kenyon et 
al. 1994), suggesting that it is much older than previously thought (see also Bertout et 
al. 1999). Mass measurements of the single GM Aur and the binary GG Tau Aa by ^^CO 
J= 2-1 line interferometric mapping have already been reported (Dutrey et al. 1998 and 
Guilloteau et al. 1999). We included DM Tau in our program because new data at 3 times 
better angular resolution would provide an independent check on the method and earlier 
results (GD98). 

Table 1 provides spectral types, luminosities, and effective temperatures for our sample. 
Most of the spectral types and stellar luminosities, L,,, in Cols. 2 and 3 of Table 1 are taken 
from Kenyon and Hartmann (1995). The luminosities are the Lj in their Table A4. All 
the luminosities in Col. 3 are evaluated at 140 pc distance. The actual distances to the 
individual stars are important for our analysis but are presently unknown; we will discuss 
this problem in Sect. 4. The effective temperatures in Col. 4, Tg//, corresponding to the 
spectral types, for all the stars, are obtained using Table Al of Kenyon and Hartmann 
(1995) for main sequence stars; we discuss the spectral-type to Tg// conversion further in 
Sect. 4. 
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2.2. Interferometer Observations 

Our IRAM interferometer observations used 5 antennas operating in snapshot mode 
(Guilloteau et al. 1992, Dutrey et al. 1996) and were carried out in winter 1997/1998. For 
BP Tau, the fuU synthesis mapping was obtained in winter 1998/1999. We observed all the 
sources simultaneously at 89.2 (HCO+ J= 1-0) GHz and 230.5 GHz (^^co J= 2-1). The 
spectra were analyzed using a correlator with one band of 10 MHz centered on the HCO+ 
J= 1-0 line, one band of 20 MHz centered on the ^^CO J= 2-1 line, 2 bands of 160 MHz 
for the 1.3 mm and 3.4 mm continuum, respectively. The spectral resolution was 0.23 and 
0.18 km/s in the narrow bands at 3.4 and 1.3mm, respectively. The phase calibrators were 
0415+379 (3C 111) and PKS 0528+134. The rms phase noise was 8° to 25° and 15° to 
50° at 3.4 mm and 1.3 mm, respectively, which introduced position errors of < 0.1". The 
seeing, estimated from observations of the cahbrators, was ~ 0.3". Baselines up to 400 m 
provided ~ 0.7" resolution for the 1.3 mm continuum data and 1.8" at 3.4 mm. The flux 
scale was checked by observation of MWC349 which has been used as calibrator at IRAM 
since 1996. Its flux is given by = 0.93 x {v/Slf-^ Jy, following Altenhoff et al. (1981), 
and subsequent measurements of the planets made at IRAM (IRAM Flux report 13). 

The GILDAS software package was used to reduce the data. At 1.3 mm, the continuum 
images are the results of summing the lower and the upper side-band data. At 3.4 mm, the 
tuning was purely in the lower side band. The continuum and line maps were produced 
using natural weighting of the visibilities. For clarity, tapering to about 1" resolution was 
used on line maps presented in Figs.l and 2, but the full angular resolution was used in the 
fitting of models by the technique. We did not subtract the continuum emission from 
the CO map because its contribution is small, but included it in the analysis. 
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3. Analysis and Results 

The essential inputs are the angularly resolved maps of the CO line emission at each 
velocity channel and the "standard" model of a rotating disk in hydrostatic equilibrium 
(e.g. DGS94). We use the procedure developed by GD98 to derive the parameters of the 
model by the method of minimization. In order to avoid nonlinear effects arising from 
deconvolution, the minimization is performed in the u,v plane. The derived parameters 
divide naturally into those relevant for the determination of the central mass and those 
that pertain mostly to the properties of the disk. Here, with the exception of DM Tau, for 
which we summarize the complete analysis, we will concentrate on the parameters required 
to determine the central mass - the disk inclination i and position angle PA^ Vo sini at To 
— 100 AU, the exponent v of the rotation curve Voir jr^'"^ and the outer disk radius Rout- 
A paper describing the disk parameters (e.g. the temperature radial distribution, the ^^CO 
abundance and spectral line turbulence width) is in preparation. 

GD98 demonstrated the parameter fitting procedure by applying it to ^^CO J= 
1-0 data for the single star DM Tau. They used D=150 pc for the distance to DM Tau. 
Since we now adopt 140 pc as the reference distance, we scale GD98's results to this value. 
On this basis, GD98's value for DM Tau's mass is (0.47 ± 0.06) x (D/140pc) Mq (Table 
3). Fig. 1 shows 7 of the velocity channel maps measured in the ^^CO J=2-l line for all 
the sources in our sample. Fig. 1 also presents velocity gradient maps which show on the 
red and blue shifted sides of the systemic velocity, the characteristic CO emission pattern 
attributable to rotation of a tilted disk. 

Fig. la shows that DM Tau, CM Aur, LkCalS and MWC480 appear to have no CO 
emission from the parent cloud. UZ Tau E and DL Tau (Fig. lb) suffer some confusion with 
the CO emission of the cloud. Confusion is maximum at the cloud systemic velocity, which 
is usually (but not always, see DL Tau) the stellar velocity. Hence it preferentially affects 
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the disk emission which is elongated along the projected disk minor axis. This results in an 
incorrect estimate of the disk inclination. In the case of CY Tau, we cannot be sure about 
the level of confusion because of the weakness of its CO emission. For UZ Tau E and DL 
Tau, the velocity channels hidden by the foreground CO cloud have been removed from the 
analysis. In the case of UZ Tau E, 4 out of 36 significant channels were removed. In the 
case of DL Tau, 5 out of 18 significant channels had to be removed because of confusion. 
This increased the uncertainty in the disk inclination quite significantly. 

The first step of the model fitting procedure is to measure the coordinates of the 
mm- wave continuum source in order to recenter the u, v plane visibilities on it. Table 2 
lists these source coordinates for DM Tau and for the other stars newly observed for this 
program. The UZ Tau E spectroscopic binary cannot be resolved by our observations; its 
continuum source coordinates must be effectively centered on the binary. 

Table 3 compares the results of the full analysis of DM Tau obtained by GD98 from 
their analysis of the ^^CO 1-0 data at ~ 3" resolution with that derived from our new ^^CO 
J=2-l data. Only the fitted parameters are shown; the gas density and ^^CO abundance 
are the same in the two analyses and do not affect the stellar mass determination because 
the ^^CO J= 1-0 and J= 2-1 lines are optically thick. We scaled GD98's values of the 
parameters that depend on physical size, and hence distance, from 150 pc distance to 
the 140 pc distance used here. The uncertainties are the formal 2a values from the 
fit. The position angles in Table 3 refer to the major axis of the disk; the PA listed in 
GD98 is that for the disk rotation axis so differs by 90° from the PAco given here. The 
agreement between GD98 and the new results is excellent for all the parameters except 
for the inclination which now has the opposite sign. It is not surprising that the higher 
resolution data better define the sense of the inclination. Its sign does not affect the mass 
determination. 
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Table 4 lists the parameters most relevant for determination of the central mass derived 
from our new observations and also from our previously published observations. Columns 2 
and 3 are described below. Columns 4 and 5 list the disk PA and i derived from the CO 
line maps. The disk outer radius, Rout, in Col. 5 is included here because its determination 
and that of i are coupled (see, for example, GD98 Sect. 4. 4). Col. 6 hsts the value of Vsini 
at the reference radius 100 AU and col. 7 the exponent of the radial power law dependence 
of the velocity. Within the uncertainties, the rotation is Keplerian. The stellar masses in 
col. 8 are derived using 



where 2.98 km/s is the circular velocity at 100 AU radius from a 1 M0 star. The derived 
masses scale as oc v'^{r)r and hence depend on distance to the star as M* oc D. 
The main source of uncertainty in the measured stellar mass is the inclination; the two 
objects in our sample with the lowest inclinations, CY Tau and BP Tau, have the highest 
uncertainties in their derived masses. This approach ignores the mass of the disk because it 
is usually only a few per cent that of the star. Even in GG Tau A, which probably has the 
most massive disk of the stars in our sample, the mass of the circumbinary is only ~ 10% 
that of the stars within it (Guilloteau et al. 1999). It is interesting that UZ Tau E has a 
circumbinary ring hke GG Tau A and UY Aur (Duvert et al. 1998), although the central 
hole is as yet unresolved. 

The 1.3 mm continuum emission is detected by a separate back-end (§2.2) and the 
detected continuum emission arises mostly from the innermost regions of the disk (see 
below) . The continuum images can therefore provide independent measurements of the disk 
orientation. Cols. 2-3 of Table 4 give the position angle, PAcont of the major axis of the 
disk, and its inclination, icont, measured from the 1.3 mm wavelength continuum image. To 
derive the inclination from the ratio of apparent major and minor axes, we used an effective 
seeing of 0.3". We did not include uncertainty in the seeing in the error estimate. The 
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agreement between the orientations and inclinations derived from 1.3 mm continuum data 
and ^^CO J=2-l data is excellent for most of the targets. The uncertainties of the values 
derived from the continuum images are larger, in most cases, than those from the CO data 
because the measured sizes of the disks in the continuum are generally smaller than in CO, 
specially when the inclinations are low. The sizes measured in the continuum are smaller 
because the dust continuum opacity is much smaller than in the CO line. The detected 
continuum emission therefore arises mostly in the dense inner regions of the disks (see, for 
example, D96). The values for CYTau suffer from low signal to noise. In the case DLTau, 
significant confusion from the molecular cloud seems to have affected the determination 
of the inclination from the CO emission (§3). The confusion affects the channels near the 
systemic velocity, which are the most important to determine the inclination. Accordingly, 
pending better measurements (e.g. ^^CO), we prefer to use the inclination determined from 
the continuum emission, ~ 45°, to derive the mass for DLTau. With this inclination, the 
stellar mass is 0.72 ± O.llM©, rather than 1.23 ± O.IIM0 when using the CO derived 
inclination. 

4. Compcirison of the Measured Mcisses and Theoretical Tracks 

4.1. The Theoretical Tracks 

The observational measure of the stellar photospheric temperature is its spectral 
type while the parameter provided by the models is usually Tg//. The surface gravities of 
the PMS stars lie between those of the main sequence stars and giants so the conversion 
between spectral type and Tg// determined for the giants and dwarfs may not apply. White 
et al. (1999) and Luhman (1999), in addition to others, have investigated this problem and 
demonstrate that the differences in temperature scales become apparent for the lowest mass 
stars, and below the hydrogen-burning limit. Luhman (1999) proposed an "intermediate" 
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temperature scale for the M stars. Its differences from Kenyon and Hartmann's (1995) 
temperature scale become significant for stars cooler than about M6. Since our sample does 
not include stars as cool as this, we use the Kenyon and Hartmann scale in the present 
work without modification. 

The HR diagram in Fig. 2a shows excerpts from PMS tracks calculated by DM97, 
BCAH, PS99 and SDF. In this and subsequent figures, the BCAH tracks for M/M© < 0.7 
are for mixing length parameter (mixing length/pressure scale height)^ 1.0 and 1.9 for 
M/VIq > 0.7 (BCAH98 and BarafFe, priv. comm.). Differences among the theoretical 
calculations for a star of given age and mass are obvious. For example, the DM97 tracks are 
hotter than those of BCAH for ~ 0.2 M©. SDF's models at 0.2 and 0.1 M© appear to 
contract more slowly than those of DM97, BCAH, and PS99. Since the different treatments 
of stellar convection, the equation of state, opacities, and stellar atmospheres become 
important in specific regions of the HR diagram, tests of the theoretical tracks require 
accurate measurements of mass over as wide a range of mass and age as possible. 

The dependence of on distance (§3) is important because while the average distance 
to the Taurus star forming region (SFR) is reasonably well known, the actual distance 
to a given star in it is not. The Taurus SFR extends for at least 15° on the sky. If its 
depth is comparable to its width, the distance to a specific member may scatter by ±20 pc 
around the 140 pc mean, a relative error of ±14%. Only CY Tau and BP Tau (Table 4) 
have a larger uncertainty of their mass measurement. For most of the stars in Table 4 the 
uncertainty in their actual mass is therefore dominated by the uncertainty in their distance. 
The distance of course also determines the star's location in the HR diagram because 

oc D'^. We will therefore present most of our results in §4 on modified HR diagrams in 
which the distance-independent parameter L/M^ is plotted versus Tg//. For reference. Fig. 
2b compares the tracks on this basis. 
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4.2. Measured Meisses on Theoretical HR Diagrams: The Single Steirs 

Fig. 3 plots the single T Tauris in our sample on modified HR diagrams using PMS 
tracks calculated by DM97, BCAH, PS99, and SDF. The uncertainties displayed along the 
horizontal axis are ±1 spectral type sub-class. Along the vertical axis, the uncertainties 
are the propagated internal uncertainties in stellar mass (Table 4) and an assumed ±10% 
uncertainty in the luminosity given in Table 1. We discuss the fit of each star to the 
theoretical tracks separately. 

Lk Ca 15: The L/M^ value for LkCa 15 hes near the 1.0 Mq tracks calculated by BCAH, 
PS99, and SDF (Fig. 3) consistent with its mass 0.97 ± 0.03 Mq at the nominal 140 pc 
distance (Table 4). These tracks yield a consistent age estimate of 3-5 MY. LkCal5's 
L/M^ values hes on the 0.8 M© track calculated by DM97. Agreement with the DM97 
track would require a distance of ~ 115 pc but would not affect the age estimate. 

DL Tau: The position of DL Tau's L/M^ in the diagrams is within one spectral type 
uncertainty of the 0.8 Mq track for each of the calculations. This is in good agreement with 
the mass derived using the continuum inclination, 0.72 ± 0.11 Mq. The theoretical tracks 
yield a consistent age estimate of ~ 2 MY. 

GM Aur: The position of the L/M^ value relative to the BCAH, PS99, and SDF tracks 
is consistent with its mass 0.84 ± 0.05 Mq at 140 pc. The tracks provide an age estimate of 
~ 3 MY. The L/M"^ value hes close to DM97's 0.6 M© track which would require a distance 
of 103 pc to force agreement with the dynamical mass. 

BP Tau: BP Tau's L/M^ values hes between DM97's 0.6 and 0.8 Mq tracks, and close 
to BCAH, PS99, and SDF's 0.8 Mq tracks, while our measured mass at 140 pc is in the 
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range 0.92 to 1.49 M©. This suggests that the distance to BP Tau may be closer than 140 
pc but not as extreme as the HIPPARCOS value, 53l}[ pc (Favata et al. 1998). The large 
uncertainty in the mass produces a large spread in the age estimate, 2-10 MY. 

DM Tau: Its dynamical mass at 140 pc, 0.55 ± 0.03 Mq, is consistent with its L/M^ 
value relative to four sets of tracks. The tracks yield a consistent age estimate ~ 5 MY. 

CY Tau: The nominal value 0.55 M© is consistent with its L/M^ relative to the four sets 
of tracks with an age in the range 2-5 MY. However, the uncertainty in its mass is so large 
that this agreement is not significant, and we cannot plot it meaningfully in Fig. 3. 

MWC480 (=HD 31648): Fig. 4 provides a similar comparisons of the dynamical mass 
of the HAeBe star MWC480. We calculated its stellar luminosity (Table 1) for a distance 
of 140 pc using the spectral energy distribution, corrected for extinction and circumstellar 
emission, derived by Malfait et al. (1998). This value, 11.5 L©, is about half that calculated 
by Mannings and Sargent (1997) from photometry available at the time. Table 1 also lists 
Grady's (1999) spectral type estimate, A4, from HST spectra, which is cooler than The et 
al's (1994) earlier value, A2-A3. BCAH PMS models do not extend to M > 1.2 M©, so we 
make the comparisons only with DM97, PS99, and SDF's tracks. MWC480's L/M^ value 
lies close to the 2.0 M© tracks for the three calculations at an age ~ 7 MY. The mass at 
140 pc distance 1.65 ± O.O7M0 suggests that MWC480 lies at a somewhat greater distance. 
A distance of 170 pc would yield a dynamical mass of 2.0 MQand would be within 2a of the 
HIPPARCOS measurement of 1311^8 pc (van den Ancker et al 1998). 
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4.3. The Binciries 

The resuhs for the binaries cannot be plotted on L/M"^ vs Teff diagrams because the 
component luminosities of UZ Tau E or component masses of GG Tau are not known. We 
discuss the binaries using conventional HR diagrams. 

UZ Tau E: Mathieu et al (1996) discovered that UZ Tau E is a single-lined spectroscopic 
binary with period 19.1 days and projected semi-major axis of the primary ai sin i* = 0.03 
AU in which is the inclination of the orbit. 

Our measurement of the total binary mass (Mi + M2) (Table 4) and the period 
determine (ai + 02) = (0.153 ± 0.003) (L>(pc)/140)V3 AU. We can solve for ai and 02, 
and hence Mi and M2 if we assume that the circumbinary disk and binary orbit are 
coplanar, ico = i*- This assumption yields Mi = (1.00 ± O.O8)(D(pc)/14O)M0 and 
M2 = (0.31 ± 0.02)(D(pc)/140) Mq. We plot the primary in the HR diagram (Fig. 5) 
assuming that the system spectral type Ml (Table 1) applies to the primary, assigning, as 
before, an uncertainty of ±1 one sub-class, and luminosity in the range 0.5 to 1 times the 
total system luminosity, 1.6 Lq. Fig. 5 shows that, with these assumptions, the primary 
lies in the ~ 0.6 region of the HR diagram for all the tracks, clearly discrepant with 
the derived mass, ~ 1 Mq. It seems very unhkely that the distance to UZ Tau E could be 
~ 60% of the nominal 140 pc. Rather, the assumption that = ico is probably wrong; 
the circumbinary disk and binary star orbits appear not to be coplanar. The spectroscopic 
binary orbit will be resolvable by the next generation space-borne astrometric missions so a 
precision measurement of the orbital parameters should be possible within a decade. 

GG Tau: White et al. (1999) compared the (L,Te//) positions of GG Tau Aa and Ab, 
the components of the ~ 0.25" separation binary at the center of the circumbinary ring we 
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observed, with PMS evolutionary tracks calculated by Swenson et al. (1994), D'Antona and 
Mazzitelli (1994), DM 97, and BCAH. They used the effective temperature scales for main 
sequence dwarfs drawn from the works of Bessel (1991), Leggett (1996), and Luhman and 
Rieke (1998). They found that the tracks calculated by BCAH were the most consistent 
with the total dynamical system we have measured and the expected coevality of the 
components. The BCAH indicated an age between 1 and 2 MY. 

Fig. 6 shows a similar comparison with the tracks of DM97, BCAH, PS99, and 
SDF. The difference between the temperature scale used by White et al and Kenyon and 
Hartmann's (1995) temperature scale is unimportant at the effective temperatures of GG 
Tau Aa and Ab. The BCAH tracks suggest primary and secondary masses of ~ 0.8 and 
~ 0.7 M0, and age 1-2 MY, in agreement with White et al's finding. The PS99 and SDF 
tracks suggest similar masses and age. The DM97 tracks also indicate a young age for the 
components, in this case less than 1 MY. The total mass indicated by DM97's tracks, ~ 0.95 
M0, is significantly lower than our measured value. To reconcile it with the measured mass 
would require that the GG Tau system lie at distance ~ 103 pc which seems unlikely. 



5. Summary and Suggestions for Future Work 

Our results indicate that: 

1) Stellar masses can be measured dynamically by the rotation of their circumstellar 
disks with high enough precision, < 5%, that meaningful tests of calculations of pre-main 
sequence evolution are possible. 

2) The BCAH, PS99, and SDF models in general are in reasonable agreement with the 
measured dynamical masses at the average distance to the Taurus SFR. To force agreement 
between the DM97 models and mass measurements for stars in the ~ 0.7 to 1 Mq range 
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(LkCal5, GMAur, and the components of GGTau) would require that these stars he at 
distances 100-115 pc which seem unacceptably near. This is a consequence of the apparently 
warmer photospheres of DM97's models for stars close to 1 Mq. 

3) Consistency of the position of the HAeBe star MWC 480 with respect to DM97, 
PS99, and SDF's tracks suggest it lies on the far side of the Taurus SFR at ~ 170pc. 

Our tests of the pre-ms tracks are limited by our present ignorance of actual distances to 
individual stars. We anticipate that this limitation will be overcome by precision distances 
that will be measured by the astrometric observatories currently under construction (e.g. 
FAME, GAIA, and SIM). All our measurements at present pertain to stars with M > O.5M0 
so we have been unable to test the models at the lowest end of the stellar mass spectrum. 
This is a serious limitation and results, at least in part, from the range of masses of stars 
formed in the Taurus SFR. We look forward to the capabihties of the ALMA mm-wave 
interferometric array which will open the rich Orion SFR to investigation. Conversion of the 
observationally derived parameter, spectral type, to that provided by the models, Te//, is a 
lingering source of uncertainty, particularly at the lowest masses. Since stars contracting to 
the main sequence represent a range of surface gravities, the conversion may depend on the 
stars age. Relief from this problem should become available soon as the model calculations 
which include stellar atmospheres provide diagnostics such as color indices and model 
spectra, for closer comparison with the observations. We look forward to a fruitful decade 
ahead. 
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Stahler, and L. Siess for helpful conversations and communicating results before publication. 
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This paper was begun during M.S.'s sabbatical visit at the Observatoire de Grenoble; M.S. 
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Table 1. Parameters of Program Stars 



Name SpT U/Lq Te//(K) Refs Refs 

(CO) (U,Te//) 



Singles: 



MWC480 


A4 


11.5 


8460 


1 


1 


LkCalS 


K5 


0.74 


4350 


2 




DL Tau 


K7 


0.68 


4060 


3 


2 


GM Aur 


K7 


0.74 


4060 


4 


3 


BP Tau 


K7 


0.93 


4060 




4 


DM Tau 


Ml 


0.25 


3720 


5 




CY Tau 


Ml 


0.47 


3720 


3 




Binaries: 












GG Tau Aa 


K7 


0.84 


4060 


6 


5 


GG Tau Ab 


M0.5 


0.71 


3800 






UZ Tau E 


Ml 


1.60 


3720 


3,7 





Note. — References for Col. 5: 1- Mannings and Sargent, 
1997, 2- Duvert et al. 2000, 3- Dutrey et al. 1996 (D96), 4- 
Dutrey et al. 1998, 5-GD98, 6-Guilloteau, Dutrey and Simon 
1999, 7-Jensen, Koerner, and Mathieu 1996 References for 
Col. 6: 1) from Malfait et al. (1998) and this work, spectral 
type from Grady (1999), 2) L* from Hartigan et al. (1995), 
3) Spectral type and L* from GuUbring et al. (1998) 4) See 
Favata et al. , 1998 and Sect. 4.1, 5) and spectral types 
from White et al. (1999). 
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Table 2. Measured mm Continuum Source Coordinates (J2000.0) 



Name 




RA 






Dec 


MWC480 




58"* 


46.264* 


29° 


50' 


37.03" 


LkCalS 


04 


39 


17.784 


22 


21 


03.46 


DL Tau 


04 


33 


39.074 


25 


20 


38.14 


DM Tau 


04 


33 


48.731 


18 


10 


10.00 


CY Tau 


04 


17 


33.727 


28 


20 


46.95 


UZ Tau E 


04 


32 


43.070 


25 


52 


31.14 


BP Tau 


04 


19 


15.835 


29 


06 


26.90 



Note. — The astrometric accuracy is better than 
0.05". 
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Table 3. Comparison of Old and New Results for DM Tau 



Parameter GD98 This Work 



Systemic velocity 


Vlsr (km/s 


6.05 ± 0.02 


6.01 ±0.01 


Orientation 


PAco (deg) 


-25 ±2 


-23 ±1 


Inclination 


i (deg) 


33 ±2 


-32 ±2 


Outer radius 


Rout (AU) 


793 ± 20 


800 ±5 


Turbulent linewidth 


A^; (km/s) 


0.08 ±0.03 


0.07 ±0.02 


Temperature law: T(r) = 


100 1, looAU ^ 




temperature at 100 AU 


Tioo (K) 


31 ±2 


32.5 ±0.5 


temperature exponent 


Q 


0.63 ±0.05 


0.63 ±0.01 


Velocity law: V{r) = Vioo 


V lOOAU / 




V sini at 100 AU 


\4oosini (km/s) 


1.11 ±0.03 


-1.17 ±0.02 


velocity exponent 


V 


0.57 ±0.04 


0.53 ±0.01 


stellar mass 


M, (Mo) 


0.47 ± 0.06 


0.55 ± 0.07 



Note. — Comparison of parameters derived from CO (1-0) (GD98) 
and CO (2-1) (this work) in DM Tau. 
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Table 4. Derived Parameters (D=140 pc) 



PAcont icont P^CO icO Rout FiooSm(i) -V 

deg deg deg deg AU km/s Mq 



Singles 



MWC 480 


-lOib 11 


26 ±7 


-32 ± 1 


+38 + 1 


545 + 5 


-2.38 + 0.02 


0.50 + 0.02 


1.65 + 0.07 


LkCal5 


62 ±5 


42 ±5 


+61 ±1 


+52 + 1 


650 + 15 


-2.30 + 0.02 


0.56 + 0.03 


0.97 + 0.03 


DL Tau^ 


44 ±3 


49 ±3 


+50 + 3 


+35 + 2 


520 + 50 


1.90 + 0.06 


0.55 + 0.03 


0.72 + 0.11 


GM Aur 


57 ± 5 


54 ±5 


+51 ± 2 


+56 + 2 


525 + 20 


2.30 + 0.08 


0.5 + 0.1 


0.84 + 0.05 


DM Tau 


-1±5 


45 ±5 


-27 + 1 


-32 + 1 


800 + 5 


-1.17 + 0.02 


0.53 + 0.01 


0.55 + 0.03 


CYTau 


-56 ±7 


47 ±8 


-30 + 7 


+30 + 10 


270 + 10 


-1.10 + 0.10 


0.50 + 0.08 


0.55 + 0.33 


BP Tau 




20 ±20 


-28 + 3 


+30 tt 


108 + 4 


1.67 + 0.06 


0.54 + 0.07 


1 24 +°-^^ 

^■■^^ -0.32 










Binaries 








GG Tau A 


7± 2 


37 ± 1 


+ 7 + 2 


+37+1 


800 


2.05 + 0.06 


0.5 + 0.1 


1.28 + 0.07 


UZ Tau E 


87 ±3 


54 ±3 


+86 + 2 


-56 + 2 


300 + 20 


-2.83 + 0.05 


0.53 + 0.03 


1.31 + 0.08 



Note. — No sign for icont ^ 1- Using the inclination derived from the continuum (see text). 
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Fig. 1. — a) Results of our interferometric observations of MWC 480, DM Tau, LkCalS, 
and GM Aur. The coordinates are arc-sec in RA and Dec centered on the peaks of the 
continuum emission (Table 2). The long and short arms of the crosses, centered on these 
positions, indicate the apparent major and minor axes of the disks. The long arms are 
scaled to a full length of 4 arc-sec and the short arms are foreshortened according to the 
inclinations. The long arms point in the direction of PA^o (Table 4). The top 7 panels for 
each object are velocity channel maps with the velocity (km/sec) with respect to the local 
standard of rest indicated in the upper left corner. The contour spacing, different for each 
object, is indicated in upper right corner of the top panel. The contour spacing is 2.7 a 
on average but varies from source to source; the smallest value is 2.2a for DL Tau in Fig. 
2. The bottom panel shows the velocity gradient map for each object, color-coded so that 
red indicates red-shift with respect to the systemic velocity. The velocity contour spacing 
(km/sec) in these panels is: MWC 480, 0.3; DM Tau, 0.2; LkCa 15, 0.3; GM Aur, 0.3. b) 
Same as Fig. la but for UZ Tau E, CY Tau, DL Tau, and BP Tau. The nearly emission-free 
channels at Yisr — 6.6 km/sec (UZ Tau) and V^s^ = 5-6.6 km/sec (DL Tau) are the result of 
confusion with the molecular cloud. The velocity contour spacing for the velocity gradient 
maps in the bottom panels are, in km/sec: UZ Tau, 0.6; CY Tau, 0.2; DL Tau, 0.4; BP Tau, 
0.3. 

Fig. 2. — Top: An HR diagram showing theoretical evolutionary tracks for stars of mass 
0.2, 0.6, and 1.0 Mq calculated by DM97, BCAH, PS99, and SDF between the ages 2 MY 
(filled circles) and 10 MY (open circles). Bottom: The same tracks as in upper panel plotted 
as L/M^ vs Te// (see text). The spectral type-effective temperature conversion is that of 
Kenyon and Hartmann (1995), see text. 

Fig. 3. — Upper Left: The results for the single T Tauri stars plotted as L/M^ vs Te// for 
theoretical evolutionary tracks for stars of mass 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 
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and 1.2 Mq calculated by DM97. The filled dots are at ages 1, 2, 3, 5, 7, and 10 MY. Upper 
Right: Same, as but for tracks calculated by BCAH. Here the indicated ages are 2, 3, 5, 7, 
and 10 MY. Lower Left: Same for tracks calculated by PS99 at masses 0.1, 0.2, 0.4, 0.6, 0.8, 
1.0, and 1.2 Mq. The indicated ages 1, 2, 3, 5, and 10 MY. Lower Right: Same for tracks 
calculated by SDF. The range of masses and indicated ages are same as for DM97. 

Fig. 4.— Top: L/M^ vs T^ff for the HAeBe star MWC480 plotted with DM97 tracks for 
stars of mass 1.0, 1.2, 1.5, 2.0 and 2.5 Mq. The dots indicate ages 1, 2, 3, 5, 7, and 10 MY. 
For the 2.5 Mq track, the isochrone dots stop at 7 MY at this age the star is on the main 
sequence. Middle: Same but for tracks calculated by PS99. The indicated ages are 1, 2, 3, 5, 
and 10 MY. The 2.5 Mq star is on the main sequence at 5 MY. Bottom: Same but for SDF 
tracks. The range of masses and indicated ages are same as for DM97. 

Fig. 5. — The primary of the spectroscopic binary UZ Tau E on an HR diagram with 
evolutionary tracks plotted in the same format as in Fig. 3. The system spectral type Ml 
is assumed to apply to the primary with an uncertainty of ±1 sub-class. Its luminosity is 
plotted in the range from total system luminosity to half this value. If the inclinations of the 

circumbinary CO disk and the stellar binary are assumed the same, the mass of the primary 
is 1.0 ± 0.08 Mq (see text). 

Fig. 6. — Upper Left: The components of GG Tau, Aa and Ab, plotted on HR diagrams 
showing the DM97 tracks. The masses and indicated ages are the same as in Fig. 3. Upper 
Right: Same but using BCAH tracks, presented in same format as corresponding panel in 
Fig. 3. Lower Left: Same but for PS99 tracks, in same format as corresponding panel Fig. 
3. Lower Right: Same, but for SDF tracks, following format of Fig. 3. 
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